Hadronic form factors and the J/ip secondary production cross section: an update 
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Improving previous calculations, we compute the D -\- D J/ip -\- tt cross section using the 
most complete effective lagrangians available. The new crucial ingredients are the form factors on 
the charm meson vertices, which are determined from QCD sum rules calculations. Some of them 
became available only very recently and the last one, needed for our present purpose, is calculated 
in this work. 



I. INTRODUCTION 



Before RHIC, RHIC physics was relatively simple. We were basically searching for a quasi- ideal gas of deconfined 
, quarks and gluons(QGP). One of the best signatures of this new state of matter was charmonium suppression [1]. 
^ ' During the last four years, due to intense work both theoretical and experimental, this naive picture changed drasti- 
— , cally. On the theoretical side, careful numerical simulations [2-5] have shown that, due to the importance of charm 
CO ■ recombination in the deconfined phase and also to final state interactions, the number of J/tp^s may stay approxi- 
mately the same. From the experimental side, especially from the analysis of elliptic fiow, came the conclusion that 
the new state of matter is not a gas, being rather a strongly interacting fluid, now called sQGP [6]. Taking the existing 
calculations seriously, it is no longer clear that an overall suppression of the number of J/ip^s will be a signature of 
QGP. A more complex pattern can emerge, with suppression in some regions of the phase space and enhancement in 
others [7,8]. Whatever the new QGP signature (involving charm) turns out to be, it is necessary to understand better 
', the mechanisms of J/^ production and dissociation by collisions with comoving hadrons. 

A great effort has been dedicated to understand J/ip dissociation in a hadronic environment. Since there is no 
O il direct experimental information on J/tp absorption cross sections by hadrons, several theoretical approaches have 
1^ ■ been proposed to estimate their values. One approach was based on charm quark-antiquark dipoles interacting with 
the gluons of a larger (hadron target) dipole. This is the Bhanot-Pcskin (BP) approach [9], which was rediscovered by 
Kharzeev and Satz [10] in the mid-nineties and updated [11,12] in the last years. Finally, the recent next to leading 
order calculations presented in [13] have conclusively shown that, for charmonium, the formalism breaks down because 
^ , this system is not heavy enough. Also considered was quark exchange between two (hadronic) bags [14,15]. The most 
' explored approach has been the meson exchange mechanism [16-20]. In our opinion the most reliable calculations of 
CTj/V'-TT done with QCD sum rules [21]. However, due to a low momentum approximation, the validity of this 
calculation was restricted to low energy reactions, close to the dissociation threshold. This is probably not enough for 
the numerical simulations mentioned above. Therefore, in order to have cross sections valid at higher energies, the 
effective lagrangrian approach still remains the most appropriate option. 

After many works on the subject some consensus has been achieved, at least in what concerns the determination 
of the order of magnitude, which, in the case of the J/iJj pion interaction, is determined to be 1 < crj/^_^ < 10 mb 
in the energy region close to the open charm production threshold. 

Once the J/V' dissociation cross section is known, using detailed balance one can attempt to estimate the charmo- 
nium formation cross section through the fusion of open charm, as, for example, DD —> J /ip -\- tt. This is known as 
secondary charmonium production. As it was first pointed out in [22], in nucleus - nucleus collisions at sufficiently 
high energies, the number of produced D and D* mesons increases and also the lifetime of the hadronic fireball 
increases. It becomes then possible that a significant number of J/^/^'s be formed by open charm fusion. Later, an 
estimate made in [23] indicated that this mechanism would be relevant only for LHC energies. The authors stressed, 
however, that their conclusion was very sensitive to the value of the J /ip formation cross section, or equivalently 
to the absorption cross section, which in that case was the one computed with the Bhanot-Peskin approach. The 
subject was left aside for some time. Recently, after the revision of the J/ip absorption cross section to larger values, 
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secondary J /ip production was incorporated in event generators in [24] and [5]. According to these simulations, the 
number of secondary J/^'s is significant already at RHIC energies. 

Given the renewed interest on the subject we shall, in this work, further refine our estimate of the J /ip interaction 
cross section, giving now emphasis to secondary charmonium production. We shall employ effective lagrangians with 
form factors calculated with QCDSR. In particular, in the present calculation we shall make use of the D*D* J lip 
form factor, which was obtained only very recently [25] and we shall also calculate the Z)*£>*7r form factor, which had 
not been calculated so far. 



II. THE EFFECTIVE LAGRANGIANS 

Since the pioneering work of MuUer and Matynian [16], there has been an intense discussion concerning the details 
and properties of the effective lagrangians which describe the interactions among charm mesons. Here we do not add 
anything new to this discussion. We shall use what we believe is the state - of - the - art lagrangian. For the sake of 
completeness and for future use we present below the effective lagrangians considered in this work: 

C-D^D. = igo^DAD^d'^TrD - Dd^'nD;) (1) 
C^D'D = 9^D'De>""^^d^'^, {daD^D + Ddc,D*p) (2) 
jC^ddt^ = igii,DDT,e^''°'^i}^duDdaTTdpD (3) 

C^DD = -ig^DD V (d^^DD - Dd^D) (4) 
Cd^D'. = -9D'D'.e^'"'^d^D:Trd^D*^ (5) 

+ {d^-iP.D*" - i^ud^Dl) D*i^ + D*!" {rd^Dt - d^^p^D*-)] (6) 
jC^D'D^ = -g^D'D^V {DnD*^ + DI-kD) (7) 

C^D'D'^ = igi,D'D'nei""^^^^,D\do,TzD*0 + ih^D'D'nei""^^d^iP,DlTTD*0 (8) 

With these lagrangians were are able to compute the processes DD — > J/'ip + tt, which involves the diagrams 
in Figure 1, the process D*D J /%p + tt, corresponding to the diagrams shown in Figure 2 and also the process 
D*D* J/tP + TT, corresponding to the diagrams in Figure 3. 

Calling pi and p2 the four-momenta of the incoming particles and ps and p4 the four-momenta of the outgoing 
particles, we can derive the Feynman rules from the above lagrangians and obtain the invariant amplitudes for each 
one of the processes in Figures 1, 2 and 3. They are given by: 



M]r = -gn-D. P3^ ( - ^ ) g^D'D e"^"'' qe Pi, (9) 



M = gD'D^P3^ — gaff X 

^ u- mfj, \ iriD' J 

g^D'De'"''^{p2-P3)eP4, (10) 
M';=g^D'D.e^p''p,,P2,P3, (11) 

Mil = 9D*DnP3. ^ 2 d'tl^DD {P2 - q)a (12) 
777-^ 
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■Mil = 9D'DnP3^ ^ 



9al3 



{P2 -P3)»{P2-P3)f3 



{Pl -P2+ P3)a- 9ul3 - {Pl + Pa) 13 9va + {P2 - P3 + Pi)u 9<jI3 



(14) 



■M^CT = -9iI>d*d-k9v 



■^lla = 9D-DnP3^ t - 9^0*0 e^"'^'' P4^ P2s 



(15) 
(16) 



D 



■^lia = -9D'D'kP3^ 2 9i)D'D e''""^ PA^Pl. 



u — m 



D 



(17) 



9tpD'D' 



t — mjj. 



QaQf} , ^ 
2^ \ 9al3 I X 



mo* 



(Q - P2)a 9ul3 + {P2 + Pi)() 9y<T - (Pl + q)y 9af3 



(18) 



M'^l, = -g^D'D' e'-^'-p,, {p2-P3)6 



u — m 



{P2-P3)a{p2-P3)l3 \ ^ 
2~ \ 9a0 I X 



D" 



mo* 



(P2 -P3- Pl)a 9^,13 + {Pl + Pi) (3 9p.o- " {P2 - P3 + Pa)^. 9^/3 



(19) 



M 



3e 



-gipD'D'iT e P3e 



n^D'D'u e PAe 



Finally the cross section for these processes is obtained with: 

da 1 1 \ - 

'dt ~ 7V 647rp2 ^ 

spin 



(20) 



(21) 



where pf a three-momentum squared in the center of mass system and the factor — comes from the average over the 

initial state polarizations. 

As extensively discussed in previous works, although the above lagrangians and amplitudes are quite satisfactory 
from the point of view of symmetry requirements, their straightforward application to the computation of cross sections 
leads to unacccptably largo results. This comes from the fact that the exchanged particles may be far off-shell and 
therefore they enter (or leave) a vertex with a very different resolving power. In one extreme case considered in the 
recent past [26], a virtual J/t/j probing a D meson, had the behavior of a parton (!). Of course, when this happens, the 
compact J/ijj almost misses the large D and as a consequence the cross section of the whole process drops significantly. 
This physics of spatial extension and resolving power is contained in the form factors. It has been realized by many 
authors that calculations with and without form factors lead to results differing by up to two orders of magnitude! 
Therefore we simply can not ignore the form factors. We must include them in order to obtain reliable results! 

Looking at the diagrams in Figures 1, 2 and 3 we notice that we need the following form factors (and the corre- 
sponding coupling constants): 



9iDUt) Cz^W 9%Ut) 9^fLM 9%Lit) giVo'it) 



(22) 



where t is the usual momentum transfer squared and in the superscript in parenthesis we denote the off-shell particle. 
This is an important distinction, because the form factors in the same vertex are very different when different particles 
are off-shell. The most reliable way to compute these factors is with the use of the QCD sum rules techniques. The 
first one of the list was calculated in [27], the second in [26], the third and fourth in [28], the fifth in [25] and they 
read: 
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si's- nit) = 17.9 ( ) = h4t,ml) (23) 



(D) 



{'S.5GeV)^-t 

Ci^W=5-8(e-(^)')=/i3W (24) 

,{t) = 20 (^e-i^y^ =hi{t) (25) 

.(i) = 13(e-(^)') =/i2(i) (26) 

5^.^.(0 = 4.8 (e(B^)) = he{t,ml) (27) 



fn}n^it)=6.2[ei^))=h,{t) (28) 
The last form factor in (22) will be calculated below. 



III. THE 7^0*0* FORM FACTOR 

In this section we shall, for the first time, compute the ttD*D* form factor using QCDSR [29,30]. In this approach, 
the short range perturbative QCD is extended by an OPE expansion of the correlators, which results in a series in 
powers of the squared momentum with Wilson coefficients. The convergence at low momentum is improved by using 
a Borel transform. The expansion involves universal quark and gluon condensates. The quark-based calculation of 
a given correlator is equated to the same correlator, calculated using hadronic degrees of freedom via a dispersion 
relation, providing sum rulc;s from which a hadronic quantity can be estimated. 

We shall use the three-point function to evaluate the D*D*tt form factor for an off-shell D* meson, following the 
procedure suggested in ref. [31] and further extended in [28]. This means that we shall calculate the correlators for a 
D* ofF-shcll and then for a n off-shell, requiring that the corresponding extrapolations to the respective poles lead to 
the same unique coupling constant. 

The three-point function associated with a D*D*-k vertex with an off-shell D* meson is given by 

^iViv^P') = I d'xd^y {Q\T{Mx)My)jl{Qm e'^'-^ e'^-^'^-y , (29) 



where = iu'fsd, = c^vU and = c-fnd are the interpolating fields for the tt", D*^, and D*~ respectively with 
u, d and c being the up, down, and charm quark fields. 

The phenomenological side of the vertex function, r^y(p,p'), is obtained by the consideration of tt and D* state 
contributions to the matrix element in Eq. (29): 

^^'^^ rnu+rnd{q^-ml.){p'-ml,){j^^-ml) 

e^^^^Vw' i + higher resonances . (30) 

To derive Eq. (30) we have made use of 

(D*-(p)]7r-(p')i5*°(9)) = i5l^22..(9')e"^^'P59Ae„(p)s(g) (31) 
where q = p — p' , and the decay constants and fu* are defined by the matrix elements 

i^\3^\^ip')) = TT^^ (32) 
m„ -I- rud 
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and 



(0|j^|£>*(p))=mc./fl.e^(p), (33) 

where e" is the polarization of the veetor meson. The contribution of liigher resonances and continuum in Eq. (30) 
will be taken into account as usual in the standard form of ref. [34], through the continuun thresholds Sq and uq, for 
the D* and tt mesons respectively. 

The QCD side, or theoretical side, of the vertex function is evaluated by performing Wilson's operator product 
expansion (OPE) of the operator in Eq. (29). Writing F^^ in terms of the invariant amplitude, 

T^4p,p') = Kp\p'\q^)e^''^'pxp',, (34) 

we can write a double dispersion relation for the invariant amplitude. A, over the virtualities and p'^ holding 
02 = _g2 £xed: 

A'-vy-.g-)^-^/;..f .. ./l;-:^^,, . (35) 

where p{s, u, Q^) equals the double discontinuity of the amplitude r{p'^,p'^, Q"^) on the cuts mg <s<oo, 0<u<oo, 
which can be evaluated using Cutkosky's rules. Finally, in order to suppress the condensates of higher dimension and 
at the same time reduce the influence of higher resonances we perform a standard Borel transform [29]: 

n(M2)^ lim hQ^r+^ (-^\\{Q^) (36) 

(Q2 = — q2 g^j^^ ^j^g squared Borel mass scale = Q^/n is kept fixed) in both variables = —p^ and 
P'^ = —p'^ M'"^ and equate the two representations described above. We get the following sum rule: 



TO„ + rud 
1 



<qq> exp(-mc/M2^ 



ds I duexp{-s/M'^)exp{-u/M"^)f{s,t,u)9{uo- u) 
1? Jo 



(37) 



where t = —Q^, 



A(s, u, t) = + u'^ + t'^ - 2su - 2st - 2tu, and Umax = s + t-ml-:^. 

We use the standard values for the numerical parameters: mu* = 2.01 GeV, = 140 MeV, = \/2 x 93 McV, 
fo' = 240 MeV, m„ +TOd = 14 MeV, ttIc = 1.3 GeV, {qq) = -(0.23)^ GeV^. For the continuum thresholds we take 
So = (tod* + A^)^ with As = 0.5 ± 0.1 GeV and mq = 1.4 ± 0.2 GeV . 

In Fig. 4 we show the perturbative (dotted line) and the quark condensate (dashed line) contributions to the form 
factor gi'^Jjj, (Q^) at = o.5 GeV^ as a function of the Borel mass at a fixed ratio M'^M^ = 0.64/(m|,. -m^). 
We see that the quark condensate contribution is bigger than the perturbative contribution for values of the Borel 
mass smaller than ^ 4.5 GeV^. However, the sum of the both contributions for the form factor, is a very stable 
result as a function of the Borel mass. The quark condensate contribution grows with while the perturbative 
contribution decreases. This imposes a limitation over the region of that we can use to study the dependence 
of the form factor. Fixing = 10 GeV^, in Fig. 5 we show, through the filled circles, the momentum dependence 

of<?(^*4.(02). 

Since the present approach can not be used at <; q, in order to extract the g-jrO'D* coupling from the form 
factor, we need to extrapole the curve to the mass of the off-shell meson D*. In order to do this we fit the QCDSR 
results with an analitycal expression. We have obtained a good fit using a exponential form: 

51^*4.(0') = 4.8 e-Q^/6 « GeV-i, (39) 

where 6.8 is in units of GeV^ This fit is also shown in Fig. 5 through the solid line. From Eq.(39) we get grrO'D* = 

dl^D'D'i.Q^ — —'m-D') = GeV~^. To check the consistency of the calculation, we also evaluate the form factor at 
the same vertex, but for an off-shell pion. In this case we have to evaluate the three-point function 
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^^;J{P,P') = J d'xd^y {0\T{Mx)My)jlm\0) e'^'- e'^^-^'^-y . (40) 
Proceeding in a similar way we obtain the following sum rule: 



ruu + irid 

- ^ r ds r du ^-s/M^^-u/M'^ ^rnMs + u-t-2ml 



(41) 



where Umtn = "z^ - and Smin = ml{l- ^g_^2 ] ■ Now we use uq = {mw + A^)^ with A„ = 0.5 ± 0.1 GeV, and 



2 rnj 
— // 

= M ^. The results are also rather stable as a function of the Borel mass. We also got a good fit of the QCDSR 
results for g^^]^.j-,.{Q'^) using a exponential form: 

5S.z..(Q') = 8-5e-«^/^-^ GeV-\ (42) 
where 3.4 is in units of GeV^. This fit is also shown in Fig. 5 through the dot-dashed line. Prom Eq.(42) we get 

g.D'D' = g'':lD'{Q^ = -^l) = 8-5 GeV-i. 

Considering the uncertainties in the continuum thresholds, and the difference in the values of the coupling extracted 
when the D* or tt mesons are off-shell, our result for the ttD*!)* coupling constant is 

5,rD*D* =8.6 ±1.0 GeV-\ (43) 

The triple vertex couplings were calculated as explained above. The quartic vertex couplings could not be obtained 
with QCDSR and we have used the prescription given in [17]: 

\/3 l\ 9aK 



^^^^-=l-6"~4J 16^ (44) 



gipD'Dn = -^gjpDD QD'Dtt (45) 

= 2 32^ ^^^^ 

where ga is obtained from: 

^'^^•^ - ^T6^ ^^^^ 

In the above expressions Nc = 3 and the triple vertex couplings are taken from our calculations. For completeness we 
present in Table I all the couplings. 





17.9 




4.0 GeV-' 


QiiDD 


5.8 




6.2 


go'D-TT 


8.6 GeV-' 


g-iliDDn 


10.0 GeV-'' 


gipD'Dn 


51.9 


giliD*D*Tr 


57.0 GeV 



TABLE I. Coupling constants used in the numerical calculations. The first four were calculated with QCDSR and the last 
three were obtained with the prescription of [17] 
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IV. THE CROSS SECTIONS 



Having all the needed form factors, we now proceed to the evaluation of the cross sections. As in previous cal- 
culations, these cross sections for secondary J/i/' production will be related to the annihilation through detailed 
balance: 

C7(3+4^1+2) - ^(1+2^3+4) (253 + 1)(254 + 1) Pj 

In figure 6 we show the J/tp secondary production cross section as a function of ^/s, without form factors. In all 
figures, the channels DD — > J/tjj + tt, DD* — > J/tp + w and D*D* — > J/ij} + tt are represented by solid, dashed and 
dotted lines respectively. In figure 7, with the help of (48) wc show the corresponding inverse reactions. As it can 
be seen, the cross sections have the same order of magnitude in both directions. Figures 8 and 9 are the analogues 
of 6 and 7 when we include the form factors in the calculations. Of course, as we stressed in the introduction, only 
these last two figures correspond to realistic numbers. The comparison of the two sets of figures is interesting only 
to estimate the effect of form factors. In previous studies doing the same kind of comparison, as for example in [17], 
the introduction of form factors reduced the cross sections by factors ranging between 20 and 50 depending on the 
channel. In that work the form factor was the same for all vertices and the cut-off, not known, was estimated to 
be between 1 and 2 GeV. Our study is much more detailed and not only each vertex has its own form factor, but, 
depending on which particle is off-shell the form factor is different. The final effect of all these peculiarities is the 
reduction of the cross sections by a factor around 7. Although significant, this reduction is smaller than previously 
expected. 

Figure 8 contains our main results. The plotted cross sections can be compared with the results of [22] and, more 
directly, with [5]. In figure 2 of [22], although the variables in the plot are different, we can observe the same trend and 
relative importance of the three channels. In that work, the results were obtained with the quark model of [15]. Our 
curves share some features with the results of [5], such as, for example, the dominance of the DD* channel and the 
falling trend of the DD* and D*D* channels. The behavior of the DD channel is quite different. In the energy range 
of \/s > 4.5 GeV our cross sections arc smaller by a factor of 2 (DD*) or 5 {D*D* and DD). These discrepancies 
are large but they are expected since in [5] all channels include the final state J/tp + p. We could not include it 
consistently because the form factors of the pDD* and pD*D* vertices have not yet been studied with our techniques 
and arc thus not available. In the model used by the Giesscn group the cross sections for D + D — > J/tp + tt and 
D + D ^ J/tp + p are similar and the same conclusion holds for the other inital state open charm mesons. If this 
would remain true in the effective lagrangian approach, then our results including both final states would come closer 
to those of [5] , giving thus a more theoretical support to the model considered there. 



V. SUMMARY AND CONCLUSIONS 



We have updated the calculations of the cross sections for J/tp dissociation and production in the effective lagrangian 

approach. The novel feature introduced in this work is the use of the form factors (23) - (28) and especially (39), 
which was obtained here from QCDSR. We believe that our results arc useful for numerical simulations of heavy ion 
collisions, such as those performed in [2-5] and also in [22-24]. The calculation of the cross sections of the processes 
considered here are complete. However, our program is not yet finished and there are still form factors to be calculated, 
such as pDD* and pD*D*. These calculations are under way and we will eventually have all hadronic form factors. 

Although no strong statement can be made without knowing the D + D ^ J/tp + p cross section, our results give 
partial support to the conclusion advanced in [5], namely that the open charm fusion cross sections are large enough 
to produce a sizeable number of "recreated" J/tp^s already in heavy ion collisions at RHIC. 
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FIG. 4. dependence of the perturbative contribution (dotted line) and the quark condensate contribution (dashed Une) 
to the gi°o'0'(Q^) at = 0.5 GeV^ The sohd hue gives the final result for the form factor. 




10 



100 



10 

aimb) 



0.1 



■ DD -> \|/+7i 
DD* -> v|/+7i 
D*D* ->v|/+7i 



4.0 4.2 4.4 4.6 4.8 5.0 5.2 5.4 

yfs{GeV) 

FIG. 6. J/tp secondary production cross section without form factors 
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FIG. 7. J/ip absorption cross section obtained through detailed balance without form factors 
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FIG. 9. J/tp absorption cross section obtained through detailed balance with form factors 
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